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Abstract 



We study expectation values of observables in three-dimensional spinfoam quantum 
• gravity coupled to Dirac fermions. We revisit the model introduced by one of the authors 

and extend it to the case of massless fermionic fields. We introduce observables, analyse 
their symmetries and the corresponding proper gauge fixing. The Berezin integral over 
the fermionic fields is performed and the fermionic observables are expanded in open paths 
and closed loops associated to pure quantum gravity observables. We obtain the vertex 
, amplitudes for gauge-invariant observables, while the expectation values of gauge-variant 

observables, such as the fermion propagator, are given by the evaluation of particular 
spin networks. 



> 

! 1 Introduction 

oo 



In [T], a method for the coupling of fermions to 3d spinfoam quantum gravity for Euchdean 
signature metrics was proposed. The resulting model was expressed as a sum over fermionic 
CO ' loops throughout the spinfoam which were then coupled to the quantum geometry via grasping 

operators [2J |3] . This yielded vertex amplitudes in the state sum given by modified 6j symbols 
incorporating the presence of matter fields. While various other methods of coupling matter to 
spinfoam models have been proposed, for example [11 [3 [6l [71 HI El [10] , the simplicial formalism 
used here, and in [11] [T^l US] , has very clear links to the classical action and closely resembles 
gauge theory on a random lattice [Tl] . 
^ ' In this paper, we firstly extend the framework of the original proposal to massless fermions. 

This will turn out to clarify, and simplify the calculation. Motivated by recent interest in the 
graviton propagator [151 [El [El [E] i we then show how to compute the fermion propagator and 
a few other simple observables. This requires a careful analysis of the symmetries of the system 
and a proper gauge fixing. Wc then integrate out the fermions using a different, more straight 
forward technique than in the original proposal. This procedure expands the observables into 
Feynman diagrams which, wc show, naturally couple to the quantum geometry. 

The organisation of the paper proceeds as follows. In section [21 we recall the proposed 
coupling of fermions. In section [3] we perform the integrations in the path integral and arrive 
at spin foam models or spin network diagrams for the various observables. In section [3] we 
demonstrate the techniques by considering some simple examples. The Appendix provides 
some useful facts on the representation theory of SU(2). 
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2 Classical theory 



Let M be a connected, oriented, compact, three-dimensional differential manifold endowed 
with an Euclidean metric g with diagonal form r/ = (+ + +). The relevant symmetry group is 
the spin group associated to the Euclidean metric r^, that is, the Lie group Spin(3) ^ SU(2). 
We will assume that the manifold M is endowed with a given spin structure, i.e. a SU(2)- 
principal bundle V mapped with a two-to-one bundle homomorphism onto the S0(3) bundle 
of g-orthonormal frames V over M. Note that the principal bundle V is necessarily trivial. 
Hence, we will choose once and for all a global trivialising section. 

Since we are focusing on fermions, we will work in the spinor representation of 3d gravity. 
To this aim, we define the linear map 7 : R'^ — >■ End (C^), mapping the Clifford algebra C(3, 0) 
(in which the (dual of the) vector space R'^ is embedded) onto the Pauli algebra, i.e., the 
endomorphism algebra of the two-dimensional complex vector space V = C^. In fact, we will 
think of the linear map 7 as taking value in the Lie algebra su(2) of traceless, anti-Hermitian 
matrices instead of the true Pauli algebra (traceless, Hermitian matrices). 

The dynamical fields of the bosonic sector of the theory arc the image e = 7(e) = e°7a, 
a = 1,2,3, of the soldering form e on V under the map 7, and a principal connection on 
the spin bundle V, that is, a spin connection uj = ^oj^CTa, where the symbols aa denote the 
Pauli matrices and the symbols -fa = ffa generate the Lie algebra su(2). The fields of the 
fcrmionic sector of the theory consist of a 2-component Dirac fcrmion tJj which is a section of 
the spinor bundle, that is, the vector bundle E ~ V x pY associated to V via the fundamental 
representation p of the spin group Spin(3) on V, and its Dirac conjugate {[> = ij:^ where the 
dagger ^ denotes Hermitian conjugation. 

2.1 Simplicial theory 

In order to regularise the functional integrals that would be used in a formal path integral 
approach to quantisation, we now introduce a cut-off on the number of degrees of freedom of 
the theory by using lattice-like methods. 

The so-called triangulation conjecture (or Hauptvermutung) states that all topological 
three-manifolds admit a piece-wise linear structure. In other words, all three-dimensional 
manifolds are triangulable. Accordingly, we can choose a fixed triangulation A of the spacetime 
manifold M and discretise the fields accordingly. To this aim, we will also need to consider 
the two skeleton k = (u, e, /) of the topological dual A* of the triangulation A, that is, the 
complex formed by the vertices u, edges e and faces / of A*. More precisely, we will work with 
the derived complex A"*" obtained by refining the dual triangulation A* by subdividing the 
dual faces / into wedges w (see [H] and references therein) . The edges of the derived complex 
A+ (which include the half-edges of A*) will be denoted e. We will choose an orientation 
together with a distinguished vertex for each wedge. The orientation of the wedges naturally 
induces an orientation of the boundary edges. This prescription is necessary to have a well 
defined model, but the end results will be completely independent of these choices. 

We finally assume that the global trivialising section of V is piecewise constant in each 
tetrahedron of A. Hence, the spinor and vector bundles associated to V are also globally 
trivialised by a piecewise constant section and we have a copy SU(2)^ of the group SU(2), and 
copies V^, V„ of the vector space V, together with its dual complex conjugate vector space 
V , above each vertex v of A*. Wc therefore have a basis (7a(w))a of su(2), a basis {eA{v))A 
of V and a basis {e^{v))A of V for each tetrahedron, or equivalently for each dual vertex v of 
A*, inside which the vector and spinor fields are globally defined. 

Simplicial fields. The fields are discretised in the standard way, namely by using their 
geometrical nature to assign values to the elements of the cellular complexes defined above. 

The (spin space) frame field e, as a one form, maps the one-simpliccs s of A in su(2). Each 
such one-simplex is topologically dual to a face / of A*, and the wedges of / are in one-to-one 
correspondence with the tetrahedra sharing the one-simplex s. Accordingly, it makes sense 
to assign the discretised co-frame to the wedges : = e'^"fa{v), where e° := e°(s) denotes 
the components of the image in of the tangent space vector s, measured with respect to 
the frame attached to the dual vertex v associated to the wedge w. Hence, for a given face 
/, the wedge variables represent the same vector measured in different frames. Accordingly, 
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these variables are related by a connection matrix in SU(2) (image in spin space of a S0(3) 
rotation) . 

This matrix g^, the holonomy P exp to of the connection along the corresponding edge e, 
is the discretised connection and is assigned to all edges of the derived complex A+ . Reversing 
the orientation of an edge maps the associated group element into its inverse : g^.-! = gl- 
Finally, the curvature around a wedge w is measured by the holonomy Gw = Yisedw 5e of the 
connection around w, by virtue of the Ambrose-Singer theorem. The matrix product is taken 
starting from the distinguished vertex and in the ordering induced by the orientation of the 
wedge. 

The fermionic fields, as sections of the spinor bundle, assign a spinor ^py and a co-spinor 
ipy to each vertex v of the dual triangulation. More precisely, to encompass the appropriate 
statistics, the components of the simplicial fermions iliyA are chosen to define a generating 

system of the Grassmann algebra 5a = A (-^ ® ^ ) associated to the triangulation A, with 

E = ®yYv and 'E* = ©^.V*. 



Simplicial action. With the fields defined above, the discretised action for the fermion/gravity 
system is the the sum of the Dirac action and an action for the gravitational field. To define 
the fermionic action, we first introduce the simplicial Dirac operator assigned to the edges of 
the dual triangulation: 

where Ue U{ge) is the holonomy of the connection along the edge e in the spinor repre- 
sentation, and Ee is the discretised version of the two-form E = e A e. It is evaluated on the 
infinitesimal triangles of A dual to the edges of A* and defined as 

= ^ ^ e^^e^^ sgn(we,u;g), (2) 

where the sum is taken over the three possible pairs (without counting the permutations) of 
wedges meeting on the edge e and on the vertex s(e), and the factor sgn(u'e, lUg) equals ±1 
depending on the sign of the associated (coordinate) area bivector. The spinor components 
are measured with respect to the spin frame associated to the source vertex s(e) of the edge e. 
The Dirac action yields 

Su[eu,,ge,'4^y,'ipv\ = ^^^Se, (3) 
or(e) e 

where the sum is taken over all orientations or(e) of all edges e of A*, and 

Se = {'lps(e),Dei)t(e)), (4) 

with .s(e) and t{e) respectively denoting label the source and target vertices of the edge e in 
the corresponding orientation, and the inner product (, ) being the standard Hermitian inner 
product on C^., i.e., introducing a basis of C^, (^,x) = ^aX^^b- 
The gravitational action is given by 

^GR [e^. , ffg] = ^ tr (e^, ) , (5) 

w 

where the trace is in the spinor representation. 

It was proven in [T] that the action for the coupled system 5gr-d ~ '5'gr + So converges 
point-wise towards the continuum Einstein-Cartan-Dirac action when the lattice spacing goes 
to zero. 

^Note this is not related to the flatness of the connection ; the vector attached to the one-simplex is in 
the orthogonal complement of the dual face and is accordingly unaffected by the rotation associated to the 
holonomy around the face. 
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Discrete symmetries. The introduction of a fixed background simplicial structure on M 
has broken the diffeomorphism invariance of the continuum theory (see [20j for a detailed 
discussion of these issues). UnUke the pure gravity case, the diffeomorphisms can not be 
recovered by a combination of gauge transformations because the introduction of fermions has 
broken the topological character of pure 3d gravity. The expectation values of observables will 
therefore be diffeomorphism covariant as opposed to invariant. The diffeomorphism invariance 
should be recovered in the continuum limit with some appropriate renormalisation as in lattice 
gauge theory. 

However, the discretisation procedure has not broken the local SU(2) symmetry. The 
arbitrariness in the choice of a basis of su(2) and of V at each dual vertex v is reflected 
in the invariance of the simplicial action under the following discrete gauge transformations 
parameterized by an SU(2)-valued sequence k : {v} — )■ SU(2): 

e^ I y ky e^ ky 

Gw k^^ Gw ky, (6) 

for the bosonic part San of the action, and 

Se I— >■ ^s(e) ^s(e) 

fPv '-^ 'tpykv, (7) 

for the fermionic sector S-d. In the first two lines, the transformation acts on the variables 
assigned to all wedges w containing the vertex v, while the transformation ([7]) takes place on 
all edges e of A*. 

Finally, note that the Dirac matrix D enjoys the following antisymmetry property: 

D^^ = ^D^A, (8) 

where the indices have been raised and lowered with the standard symplectic metric eAB on 
C^. Our conventions arc the following. The two-dimensional totally antisymmetric tensor is 
normalised such that eoi = e*'^ ~ +1, which implies that eab^^^ ~ ~^a- The raising and 
lowering of indices occurs according to 

= e^'^Xs, and U = ^""^BA = -^ab^"" ■ (9) 
This implies in particular the following property for multi-component spinors 

m-^:::a... = -m,,,^:::^-. (lo) 

Simplicial observables. In this paper, we are interested in the study of simplicial fermionicH 
observables O/.a := Of, that is, polynomial functions of the simplicial fermion fields. In 
particular, we will study the following observables 

• C'y"^(a;, y) — tpys '■ the gauge- variant two-point function of the fermionic field. 

• Of{e) = (V'sCe)) '^e''Pt(e))-i with e the path linking s(e) to t{e) and T4 the corresponding 
holonomy. We will call this gauge-invariant observable the Polyakov line. Note that 
this observable is sometimes argued to be the correct physical observable for the fermion 
propagator (see [H] and references therein). 

• Of{e,e') = M{e) M{e'), with M{e) = {ips(e), ^e'0t(e)); the two-point function of the 
Polyakov line. 

With the discretisation outlined above, the naive definition of the expectation value of an 
arbitrary simplicial fermionic observable Of yields 

(0/)oa-o = ( n / dej [U [ dgA ( [ dfii^y, ijy)) Of e'^«--. (11) 



^These observables are called fermionic, even if they obey bosonic statistics, because tliey are functionals 
of the spinor fields. 
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Here de^ is the Lebesgue measure on su(2), dg^ is the normalised Haar measure on SU(2), 
and the symbol c?/i('!/'„, tpv) denotes the Berezin integral on Q/\. The latter measure is the 
element of Q^^ defined by 

1 1 

dniip^, ipv)= n n 

V—n V—n 

where n is the number of vertices of A* and 
Finally, the normalisation faetor 



n / f n / i I 



2oH-o =(11/ deu, 1(11/ dg,]( I V-.) ) e'^-^-- (12) 



is the diseretised path integral for the coupled system. 

Gauge fixing. Aecording to a theorem in lattiee gauge theory by Elitzur [52], the expecta- 
tion value of a locally gauge dependent observable, such as the fermion two-point function, is 
zero without gauge fixing. This result is reproduced in the spinfoam formalism as the expecta- 
tion value of a non-gauge fixed, locally gauge-dependent observable will result in spin network 
vertex amplitudes with a single open end which vanishes by Schur's Lemma (see equation ([M)) 
in the Appendix). It is therefore crucial to appropriately fix the above gauge freedom in the 
naive expression ([TT|) . This can be achieved following the next lines [231 \E\ . 

The idea is to introduce a maximal tree of the triangulation's dual T C A*, that is, a 
sub-complex of the one-skeleton of A* touching every vertex of A* without ever forming a 
loop. One can pick an arbitrary preferred vertex r on T, called the root, and denote u the 
endpoints of the tree (recall that T is an open graph). We will note At the simplicial pair 
(A,T). Using this tree and the gauge invariance, it is possible to reach the gauge = 1 
for all e in T. Remarkably, the Faddeev-Popov determinant associated to this gauge-fixing 
procedure : 

^FP=]l[ n ^iKL 9e ktie}), (13) 

with dk denoting the normalised Haar measure on SU(2), is equal to one [211 H] because 

- 1, (14) 

where 7(u) is the unique path linking the endpoint u to the root of the tree. 

If the observables are gauge-invariant, one simply follows the standard Faddeev-Popov 
procedure associated to the above gauge-fixing prescription. If the observables are not gauge- 
invariant, the situation is more involved. One procedure to compute a gauge- variant observ- 
ables in lattice gauge theory involves averaging over the gauge group to define a gauge invariant 
function which can then be computed |25| . Alternatively, following Faddeev and Popov, one 
can take the definition of a gauge dependent operator to include the necessary gauge fixing 
[26] 

(^^/>oB.-o = ^ ( n / dej [II [ dgA ( [ d^i{i^,, V.)) n ^^9e) Of e'^^--, 

(15) 

where the normalisation factor is the path integral (|12p gauge fixed by the insertion of the 
simplicial gauge fixing function nggT<5(ge)- For a gauge invariant operator, this reduces to 
equation (jlip . while it defines the expectation value for a gauge- variant observable. One can 
show that the expectation values of gauge invariant observables are independent of the choice 
of gauge fixing tree T. The gauge variant observables necessarily depend on the choice of tree. 

This describes the classical simplicial theory. We can now proceed to the quantum com- 
putations. Indeed, having reduced the number of degrees of freedom from infinite to finite by 
the introduction of a discretisation scheme, we can now compute the integrals and accordingly 
access the quantum aspect of the simplicial theory. 
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3 Quantum theory 



To compute the integrals in (fTSl). the idea is to first integrate out the fermionic degrees of 
freedom, in other words, perform a Feynman diagram expansion of the matter sector, before 
computing the gravitational integrals. As a result, (C'/)gr-d is expressed as the expectation 
value of a bosonic observable Ob of the pure quantum gravity theory : 



Therefore, the expectation value of a fermionic observable can be written 



where 



and, using the fact that the Dirac action Su is an element of ( 



(^J dfiii:^, e'^° =detD. 



The next step is to expand the bosonic observables Ob in fermionic paths. 



(16) 



Ob^[ I d^ii^,, V.) ) Of e^^°, (18) 
'Sa 



3.1 Feynman diagram expansion 

The advantage of first performing the matter integrals is that the Feynman diagrams define, 
as mentioned above, pure quantum gravity observables that naturally couple to the dynamical 
geometry. The essence of the expansion relies on the properties of the Berezin integral. The 
techniques used here are different than the ones introduced in [T], which relied on the intro- 
duction of symplectic Majorana fermions. Here, we explicitly integrate over the Dirac fields 
which offers a more direct way to compute the integrals. 

The first step consists in Taylor expanding the exponential of the Dirac action. We use the 
fact that each term in the sum ^ is an even element of the Grassmann algebra, that is, has 
overall bosonic statistics SeSe' = Se'Se, for all couples of edges e, e', and that the expansion is 
finite because (Se)^ = 0, for all edges e and all N > 2. This leads to the following expansion 
of (Ull) 

^'^{ig '^f'^'^"' '^^'^) '^Z ri(l+^«^e- (^l + ja5e-i-y(5,-0'), (19) 

where a ~ 8(27rG)^ (we have rescaled the triad such that e i-^ e/lQnG). The above expansion 
produces a finite but possibly large number of contributions. However, the number of terms 
surviving the Berezin integral is drastically reduced. 

Being a derivation, the Grassmann integral yields a non-zero result only if the integrand 
is the top form on Q/^, or in other words, if for each vertex v of A*, there is a product of two 
non-identical spinor and two non-identical co-spinor components. 

Consider first the vertices where there are no field insertions. The edge action Se appearing 
in the expansion has a co-spinor (resp. a spinor) siting at the source (resp. target) of the edge 
e. Hence, the only contributions surviving the integration are such that, for each such vertex 
V, there is a product of four edge actions Se such that v is exactly a source vertex for two edges 
and a target vertex for the two other edges. Note that the four edges need not be distinct. 

To integrate such contributions, we firstly consider the Berezin integral for the Grassmann 
algebra associated to a single vertex. Dropping the vertex label for a moment, it is easy to 
prove the following identity 

J d/i(V^, ^) i^"" = -e^'^'eAc- (20) 
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If one implements this formula at a vertex v satisfying the requirements discussed above, it 
is immediate to see that the integration connects, via the symplectic metric, the indices of 
the Dirac matrices De (appearing in the edge actions Se) associated to the edges meeting at 

V pairwise : the two matrices associated to the two ingoing edges are paired together and 
likewise for the two matrices corresponding to the two outgoing edges. This implies that two 
matrices being connected are necessarily associated to edges having opposite orientations. The 
symplectic metric raises (resp. lowers) the second (resp. first) index of one of the two ingoing 
(resp. outgoing) matrices. This index then gets contracted with the first (resp. second) index 
of the other ingoing (resp. outgoing) matrix. The orientation of the edge associated to the 
matrix with raised (resp. lowered) indices is flipped using the symmetry ^ for the orientation 
of the edges to be consistent with the matrix product. Therefore, the Berezin integral at vertex 

V yields two products of Dirac matrices (one for the ingoing and one for the outgoing edges) 
with consistent orientation of the associated edges. Repeating the procedure for each vertex 
with no field insertion will accordingly produce products of closed sequences of Dirac matrices 
associated to edges forming loops in the dual triangulation A*. The number of Dirac matrices 
appearing in each such sequence will necessarily be even because the matrices getting paired 
are associated to edges with opposite orientation. 

Let now w be a vertex with a fermionic field insertion. If a spinor (resp. a co-spinor) 
component is inserted on a vertex v, one only requires another spinor (resp. co-spinor) com- 
ponent and two co-spinor (resp. spinor) components to saturate the vertex. Therefore, the 
contributions surviving the integral are the following. For each vertex v where a spinor (resp. 
a co-spinor) is inserted, there is a product of three edge actions Se where the corresponding 
edges, which again are not necessarily distinct, are such that v is exactly a target for one edge 
(resp. two edges) and a source for the other two edges (resp. for the other edge). This implies 
that for each spinor (resp. co-spinor) insertion, there must be also a co-spinor (resp. spinor) 
insertion on some other vertex to obtain a non- vanishing result. 

The integration at a vertex v where a spinor (resp. co-spinor) is inserted then connects 
the two Dirac matrices associated to the two outgoing (resp. ingoing) edges as before, but 
the Dirac matrix corresponding to the ingoing (resp. outgoing) edge is left open with a free 
second (resp. first) index corresponding to the index assigned to the spinor (resp. co-spinor) 
insertion. The other index of the open Dirac matrix gets connected at one of the adjacent 
vertices with another matrix and, using the orientation flip ([S]), we obtain a product of Dirac 
matrices along an open oriented sequence of edges which starts (resp. terminates) at v and 
terminates (resp. starts) where another field insertion given by a co-spinor (resp. spinor) has 
occurred. This path is necessarily supported by an odd number of edges. 

Therefore, the result of the Berezin integral with p spinor and p co-spinor insertions on the 
vertices xi, X2p of A* can be stated as follows. The following applies for observables with at 
most a single field insertion on each vertex, but can be trivially extended to higher insertions. 
Each non-vanishing contribution is labelled by a graph r(j;i, ...,X2p). This graph is the union 
of 1) a union of oriented loops C with 2) a union oi p oriented open paths V connecting the p 
vertices with spinor insertions to the p vertices with co-spinor insertions. 

This graph is such that 

1. It has overall lengtf|f| Er given hy 2n — p 

2. The length Ec of each connected loop component C is even, while the length E-p of each 
path components V is odd 

3. Each vertex v of the dual triangulation A* is either traversed twice by the edges of F, 
either traversed once if v admits a field insertion, i.e., if w € {a;i, ...,X2p} 

4. Each open path V is oriented such that 5(7^) (resp. t{V)) is a vertex with a spinor (resp. 
co-spinor) insertion. 

5. Consecutive edges of the graph F cannot go back and forth along the same dual edge 
e € A* unless that particular connected component of F is of length two. 

A graph satisfying all these properties will be called admissible. 

All possible, admissible graphs can be obtained by the following graphical method. Draw 
all possible sets of arrows of cardinality 2n connecting all the vertices of the dual triangulation 

^The length of a graph is the number of edges of the dual triangulation supporting the graph. 
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such that each vertex without field insertion has exactly two ingoing and two outgoing arrows, 
while a vertex with a spinor (resp. co-spinor) insertion has two outgoing (resp. ingoing) and 
one ingoing (resp. outgoing) arrow. For each such set of arrows, the corresponding admissible 
graph is obtained by connecting at each vertex the ingoing arrows with the ingoing arrows and 
likewise for the outgoing arrows. The vertices with field insertions will accordingly be sources 
or targets for open lines. The connection process flips the orientation of one of the two arrows 
being paired to obtain a consistent orientation. For example, on a simple four point lattice 
with a field insertion 'ipi'ip2y the connecting of the arrows occurs as follows 



1 



1 





A 



It may seem that there is an ambiguity when choosing which arrow to flip; the two possible 
choices correspond to the two possible orientations for the loops and paths. However, the two 
contributions corresponding to the two possible orientations of a loop are in fact equal because 
of the orientation symmetry ([8|). For a path, the boundary data flxes the orientation such that 
the path always starts at a spinor insertion and terminates at a co-spinor insertion. Thus, 
there is no ambiguity. 

We are now ready to establish a set of 'Feynman rules' to obtain the bosonic observables 
Ob in terms of the fermionic ones Of. These rules are as follows. 

• Consider all possible admissible graphs determined hy Of . 

• To each admissible graph T{xi, X2p), assign respectively to each open path V and to 
each closed loop C of r(a::i, X2p) the ordered products of Dirac matrices Dg 

Dr = (m)^^ Y[ De, Dc = (m)^^ tr [] D^, 

following the orientation of the path V or the loop C to which they are associated. 

• Assign an global sign £t{xi,...,x2p) = to each graph r(a;i, a;2p) which controls the 
overall sign appearing as a result of raising and lowering some of the indices with the 
symplectic metric (see equation ([5])), using the orientation antisymmetry ([H]), the 'see 
saw' property pop . and rearranging the fermionic fields to compute the full Berezin 
integral. 

• Sum over all such contributions: 

= X! ^rixi.,....x2p)Irixu...,x2p), (21) 

r(a:i,...,a;2p) 

where Ir{xi,...,x2p) is the amplitude associated to the admissible graph r(a;i, ...,X2p) con- 
structed from the above assignments, and, depending on the obervables, the holonomies 
Ve appearing in the definition of O/ . 

One may wander if there are no numerical factors supplementing the signs in the rules 
depicted above. The only factors to take into consideration are the 1/2 factors coming from 
the quadratic terms in the Taylor expansion of the exponential. In fact, these factors cancel 
exactly the multiplicities of the contributions to which they are associated. Indeed, the 1/2 
factors are associated to loops going back and forth along a single edge because they correspond 
to terms of the form (S'e)^ in the Taylor expansion. But for each such term there is also a 
contribution of the form (S'g-i)^ giving rise to exactly the same loop contribution which hence 
appears twice. 

Using this framework, we can study the expectation value of all polynomial observables in 
the fcrmion fields by computing the Berezin integrals according to the rules depicted above. 



8 



Firstly, wc can consider the somewhat trivial ease where there are no field insertions p — 0^ 
i.e., the computation of the functional determinant detD. In this case, the open paths P are 
simply reduced to the null graph and the admissible graphs are accordingly given by collections 
of closed loops which correspond to the terms of order 2n in the expansion parameter a. An 
example is given by a single loop going through all the vertices of the triangulation twice. 
Another example is provided by a chain-type graph given by a collection of length two loops 
each going back and forth between two adjacent vertices. The amplitude Ir associated to a 
such admissible graph F is thus given by 

Ir^llDc, (22) 
c 

where the product is over all the loops £ C F. We therefore obtain that the functional 
determinant yields 

dctD = Y,^rl[Dc. (23) 

r £cr 

We can also compute the bosonic observable associated to the gauge-variant two-point 
function of the fermionic field 

The corresponding bosonic observable is defined by specifying the set of admissible graphs 
F(a;i,a;2) and the corresponding amplitudes It{xi,x2)- The admissible graphs are characterised 
by the fact that they contain one open path V{xi,X2) going from xi to X2- An example is 
given by an open path and a single closed loop satisfying the admissibility requirements. The 
amplitude associated to an admissible graph F(a;i,a;2) yields 

Irix.^x.fB - D^ix,,x:)B {I[Dc^ . (24) 

Note that, for all admissible graphs F, the above amplitude is of order 2n — 1 in the expansion 
parameter a = i/rn^, where ra^ = 1/2t:G is the Planck mass. Therefore, since the two-point 
function is given by the quotient (|17p of the expectation value of the above bosonic observable 
with the expectation value of the functional determinant, it will be of mass dimension two as 
expected for the three-dimensional fermion propagator in the spatial representation. 
The bosonic observable corresponding to the gauge-invariant Polyakov line 

Of{e) = (-03(e), Vetj}t(e)), 

is associated to the same type of admissible graphs as the two-point function because the field 
insertions are of the same type. The corresponding amplitudes arc given by the expression 



(25) 



where V{x2,xi) is a path connecting the endpoints X2 = t(e) to xi = s(e) of the Polyakov 
line. 

Another interesting observable is the two point function of the Polyakov line 

C/(e,e') = (V's(e), Ve1ptie)){llJs{e'), K' tptie'))- 

If xi, X2 and ^3, X4 respectively label the source and target vertices of e and e', the result of the 
integration over the fermion fields is labelled by two types of admissible graphs F(a::i, ...,2:4). 
The first type contains two open paths linking X2 to xi and 2:4 to x^. The corresponding 
amplitude yields 

= tr ye^7'(a;2,a;i) tr Ve'Z3p(2:^_2;3) ^ ^£ V (26) 
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The second type of admissible graphs admit two open paths respectively linking X2 to and 
2:4 to xi, and the corresponding amplitude follows 



h{x^,...,xi) 



The above expression suggest that the amplitudes associated to both types of gauge in- 
variant observables considered here can be written generically in the following form 

•PG7 / ^ ' 

Here, the first product is over all possible closed loops 7 formed by composing the open paths 
V of the graph F with the edges e supporting the holonomies 14 defining the observable O / . 
The edge 6(7-") denotes the edge e supporting the observable holonomy which is such that 
t{e)=B{V). _ 

The obtained amplitudes are not yet expressed as products of holonomies as one would 
expect. In fact, all the above amplitudes are given by a sum over 2^^ contributions since 
each Dirac matrix is a sum of two terms (see equation ([T])). All these contributions are 
associated to the same graph F but vary by the pattern of the assignments of holonomies 
and simplical area two-forms to each edge of F. We will call each such term a configuration 
associated to the graph F and note it /p, with the parameter c = c(F), ranging from to 
2^"" — 1, labeling the configuration. An example of a such term can be obtained for each 
observable given above by replacing the Dirac matrix Z?e with the quantity Sgt^e- Another 
configuration is obtained, for example, by assigning the factor J/gSe-i, instead of Set/e, to one 
of the edges. Hence, for each graph F, /r = -^r ^"^^ ^ bosonic observables will be written 
as a sum 

^f = E^r/r- (29) 

r,c 

We close this section on the fermionic integration with a remark. Interestingly, considering 
the fermionic observable given by the top form 



leads to a bosonic obervable simply given by 

a = ±l. 

This implies that the expectation value of this observable is simply given by the ratio of the 
Ponzano-Regge amplitude and the expectation value of the fermionic determinant. We now 
turn towards the integration over the quantum geometry fluctuations. 

3.2 Gravitational integrals 

We have shown that the bosonic observable Ob corresponding to an arbitrary fermionic ob- 
servable Of could be written as the finite sum (f29l) . where /p = I^(g.^,U^ As a result, the 
expectation value of a fermionic observable O / can be written as a finite sum over spinfoams : 

^Of)..-. = ^^Y^^^.A^. (30) 
where the amplitude Ap associated to the graph F and to the configuration c is given by 




The calculation of the observable now boils down to evaluating each term pip appearing 
in the sum over spinfoams (|30l) in a similar way to evaluating Feynman diagrams in ordinary 
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quantum field theory. To this aim, we now need to perform the integration over the discretised 
frame field and holonomy variables. 

Let us take a closer look at the functional I^{eu,,Ue) associated to a particular term F in 
the Feynman diagram expansion and to a particular configuration c. The amplitudes /p for all 
observables considered here share the same structural properties : the functional /p(e„,, Ue) is 
given by products of terms of the form SgC^e or Ue'S^-i along the edges of the graph F which 
can close or not. There can also be holonomy matrices Ve linking the open ends of the graph. 
We therefore need to understand how to integrate the Dirac matrices and holonomy variables 
against the complex exponent of the gravity action. We start by considering the case of gauge 
invariant observables (^5]) . 



3.2.1 Gauge invariant observables 

To calculate gauge invariant observables, we will use the fact that, because the Fadeev-Popov 
determinant is equal to one, the gauge group is compact and the Haar measures used in the 
path integral are normalised, the expectation value of a gauge invariant observable calculated 
by inserting the gauge fixing function JleeT ^i9e) in the path integral is equal to the expectation 
value obtained without the gauge fixing. Therefore, for the computation of gauge invariant 
observables, we will simply omit the gauge fixing function. 

To understand how to perform the gravity integrals, we will choose a particular configura- 
tion /p associated to an arbitrary graph F in the sum over configurations appearing in ()28p . 
It will be sufficient to consider only one type of configurations as the techniques developed for 
this particular case will equally apply to all other terms. This particular prototypical contri- 
bution is the one obtained by assigning only terms of the form SgC^ to the edges of the graph 
F. We arrange the labels c such that this configuration corresponds to c = 0. It is then given 
by the expression 

Integration over simplicial triads. The Dirac matrices ((!]) involve holonomies and dis- 
cretised area two-forms ^ which are quadratic in the simplicial triad field. To integrate 
over the wedge variables e^,, the idea, borrowed from conventional quantum field theory, is 
to replace polynomial combinations of the discretised triad appearing in the path integral by 
source derivatives [2 13] 

i S 
^ 2 5 Jul ' 

acting on an appropriate generating functional [1], where the sources J„, = J^,7a represent 
a Lie algebra valued two-form J discretised on the wedges w. The role of the free theory 
is played by the pure topological 3d gravity action. In this framework, the discretised area 
two-form Ee = "^ela appearing in the Dirac matrices becomes the differential operator 



= E ^"'''A:|rZ^«gnK,u;^), (33) 

and the amplitude pip is re-writen as 



= (ia) 



7 ve-yeev ] \ c eec / 



(34) 

,7=0 

where the generating functional, after integration over the simplicial triad, is given by 
\ g "'su(2) J J y ^ J ^ 

(35) 

where (a, b) is a pair of cumulative indices such that a runs over all of the in each path V 
in 7, and b runs over each b^ associated to each edge e of each loop configuration C. The delta 
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function, which appears because the action appearing in the generating functional is hnear in 
the simphcial triad, is assumec0 to be on the Lie group SU(2). 

The next step consists of expanding the delta functions in terms of characters. Let j £ N/2 
label the unitary, irreducible representations : SU(2) Aut Vj of the Lie group SU(2). The 
Peter- Weyl theorem states that 

V5eSU(2), = ^ dimjx,(5), (36) 

jGN/2 

where dim j = 2j + 1 is the dimension of the representation j and Xj = tr tt^ is the correspond- 
ing character. This step transfers the calculations to the arena of the representation theory 
of SU(2). The idea is then to expand the characters appearing in the generating functional 
(|35p into products of matrix elements of the group variables associated to the boundary of 
the wedges. The same group element then appears into a number of different representation 
matrices and implementing the integrals over the simplicial connections then boils down to 
understanding how to integrate the tensor product of representations. All the formulae given 
in the following section are proved in the Appendix. 



Group integrals. The group integrals are easier to understand with the aid of the graph- 
ical calculus of spin networks. The graphical methods offer the advantage of accurately per- 
forming representation theory calculations which would involve lengthy formulae if written 
algebraically. 

A spin j representation matrix tt-' evaluated on a group element g is depicted as an oriented 
line coloured by the spin j traversing a little circle corresponding to the group element g : 

j 



The trace of a such matrix is simply obtained by joining the two ends of the open line to form 
a closed loop. If the group element is the identity, the circle is simply omitted which implies 
that a spin j line without a circle is the identity operator on Vj. Each line comes with an 
arrow, which is inherited from the orientation of the derived complex A"*", and reversing the 
arrow is equivalent to switching from a representation to its dual. Here it is assumed that, 
when no arrows are explicitly given, the orientation flows from the top to the bottom of the 
diagram. 

The Haar integral over gi of a representation matrix tt-' (g) is understood when the circle is 
replaced by a rectangular box (a cable). Thus, the Haar integral of the tensor product of n 
representation matrices with spins ji, j„ evaluated on the same group element 



dg TT^^ {g)® TT-'" (5), 

/SU(2) 

is represented by a box in which the corresponding n lines enter and exit 

h ■■■ jr. 



(37) 



(38) 



A box defines an endoniorphism of the tensor product of the corresponding representations 
which projects onto the trivial subrepresentation. Accordingly, the result of the integration 
can be expressed in terms of intertwining operators. The diagrammatic representation of an 
intertwining operator t of unit norm between the tensor product of n representations ji, jn 
and the complex numbers 

i e Hom(V,-i (g) ... (g) V,„,C), (39) 

*In fact, the naive procedure followed here yields [5] the delta function on SO(3) and not the one on SU(2). 
To obtain the delta function on SU(2), one needs to insert the observable Yl^ 5 (1+ ^tr Gw) in the path integral 
[27) . In what follows, we will assume that this trick has been implemented and that we are thus working with 
a delta function on the spin group. 
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is given by a vertex, or node (to avoid confusions with the vertices of the dual of the triangu- 
lation), with n ingoing hues: 



Ji \/ Jn (40) 

+ 

For this diagram to accurately describe an intertwining operator, one needs to specify an 
ordering of the lines to match the order of the tensor product in the corresponding algebraic 
expression. Our conventions are such that the order of the tensor product is read clock-wise 
(resp. anti-clock-wise) starting from the line marked with a bullet point if the vertex carries a 
+ sign (resp. a — sign). Therefore, the picture (|40[) is precisely the diagrammatic form of the 
intertwiner t in p9|) . The reversal of one or more arrows in the above diagram would represent 
an intertwining operator between the tensor product of the representation associated to the 
ingoing lines to the tensor product of the representation assigned to the outgoing lines. 

We are now ready to perform the group integrals. There are three different types of Haar 
integrals appearing in the amplitude depending on what type of edge the integration is assigned 
to. The first type is when the edge e belongs to the interior of a face / in which case it is shared 
by precisely two wedges and the corresponding group element thus appears in two characters. 
For every such interior edge e in a given face /, we use the orthogonality of the characters 

S 

dg X]{gig)xk{g^^g2) = :tt^Xj(.9i52)- (41) 

su(2) dimj 
Diagrammatically this is 
j 





(42) 



Implementing the above formula repeatedly forces all the spins associated to wedges be- 
longing to the same face / to be equal — jf for all w in f. The various powers of the 
dimension factors dimj = 2j + 1 appearing when performing the integrals combine with the 
dimension associated to each wedge by virtue of p6p and produce, for each face /, a weight 
given by {2jf + l)^f where Xf is the topological invariant Euler characteristic of the face /. 
From here on, we will assume that each face is homeomorphic to a disk and thus that x/ = 1 
for all /. The generating functional ([55]) thus reduces to 

Ma,b) / n - TT ^™ I 



with 



3f f 



in I dgA n n i-.u,y<v) f n*^- n ^^^-^a (44) 

V e •^su(2) ) y ^ ve^eev J \c eec J 



AT'''\jf,J) 



X 

/ 

Here, the source insertions JJ are along all the vertices w of a given face / (m := uif denotes 
the number of vertices, i.e the number of wedges contained in /). The edges e of A+ that we 
have not integrated out, that is, the edges belonging to the boundaries of the faces / (which 
are half-edges of A*), have been recombined into edges e of A*. 

To calculate (1441) . we need to perform the remaining integrals which, as we have seen, 
are associated to edges of the dual triangulation A*. The edges of the dual triangulation 
are shared by three dual faces. If a such edge e does not support an edge of the graph 
r, the corresponding group element appears in three characters and the associated integral 
involves the tensor product of three representations. For every such edge we use repeatedly 
the following formula 

/ dg tt' (g) ®TT^{g)®Ti^{g)^ ii\ (45) 

JSU(2) 
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where t : Vi ® V,- ® 14 — > C is the Wigncr 3j map, while t* : C Vi <Si Vj <Si Vk denotes the 
dual map. Introducing a basis {e^a)a=-j....j of Vj, the following evaluation 



i j k 

a /3 7 



(46) 



is the standard Wigner 3j symbol (see the Appendix for an account of our conventions). In 
the diagrammatic language, the formula (05]) reads 

i j / k 



J 




(47) 



Note that no marking of lines is required for three- valent intertwiners, we only need to specify 
an orientation of the vertex. This is because of the cyclic symmetry of the 3j symbols (see 
equation in the Appendix). 

Note that we have only given the integration formula for the case where all three lines are 
pointing in the same direction, from top to bottom. If one line coloured by, say spin fc, is 
pointing from bottom to top in equation (j87p . the formula must be corrected by multiplying 
the right hand side by a phase given by (— 1)^*^. This factor is the norm (squared) of the 
corresponding intertwiner. Since the Haar measure in invariant under the transformation 
g we now have covered all possible orientations configurations; if a line is pointing in 

a direction opposite to the two others, we pick up a phase given by —1 to the power twice the 
spin colouring the line. 

Now, if an edge e does support an edge of the graph F, the corresponding integral involves 
the tensor product of the three representations assigned to the three wedges meeting on e but 
also one (resp. two) fundamental j = 1/2 representation associated to the holonomy matrix 
(matrices) Ue appearing in the Dirac matrix (matrices) attached to the edge e. If there is 
only one spinor representation associated to e, the relevant formula to use is that of the tensor 
product of four representation matrices 



SU(2 



dg U{g) ® 7T\g) (g) (g) ® n'^ig) =^is C 



(48) 



where the spin s labels a basis of the vector space of 4- valent intertwiners obtained by choosing 
a recoupling scheme and decomposing the tensor product of two out of the four representations 
using the Clebsch-Gordan intertwining maps (see Appendix). For instance, if we choose the 
recoupling 

Hom(14 «) V,- «) Vfc Vi, C) = (Hom(T/, ® Vj, Vs) (g) Rom{Vk ® Vi, V,*) 



for arbitrary representations, the four-valcnt intertwiner is related to the 3j symbols by the 
following evaluation 



i j k I 

a l3 J S 



voir 



« J 
a (3 



s k I 
e 7 (5 



(49) 



where the indices are raised and lowered, i.e., the lines inside a given column are interchanged, 
according to the convention displayed in the Appendix and summation over repeated indices is 
understood. One could obviously choose another recoupling of representations. The coefficient 
appearing in the change of basis between the chosen decomposition and the one given by a 
different recoupling of representations is (proportional to) a six-j symbol (see equation (jlOOp 
in the Appendix). 

In the diagrammatic language, equation (j48|) yields 



1 2 

t 

1 


t 


J 




1 
1 
1 







= E 



.1 

2 



J, 




+ S 

_ s 



= dim i 




(50) 
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Note our convention in which the fundamental representations are represented by dashed hues. 
As in the three- valent case, the formula would need to be corrected by a factor (—1)^' if the 
line of spin I was pointing in a direction opposite to the three others, and by a factor (— 1)2('=+0 
if the lines k and I where going in opposite directions with respect to the two others. 

There can also be two fundamental j = 1/2 hues El associated to the edge e. The corre- 
sponding integral involves the tensor product of five representations and the relevant formula 
is 

dgUig)®Uig)®T:'{g)®T:^g)®T:''ig) = Y,^^,t ^-It: (51) 



SU(2) 



where the spins s, t now label two channels in a particular recoupling scheme. We will choose to 
couple the two spinor representations together and likewise for the j and k representations. The 
first coupling gives rise to s while the second produces the label t. For general representations, 
but under the same recoupling scheme, the five-valent intertwiner L^.t is also related to the 3j 
symbol as follows 



i j k I m 
13 S C V K 



= Vdi 



i dim t 



s,t 



« 3 M 
/3 6 s 



ska 



t I m 
a T] K 



(52) 



Therefore, the diagrammatic representation of ([5T|) yields 



1 1 

1 2 1 2 


i 


3 


k 


+ + 
1 1 










= dim s dim t 


1 1 
1 1 
1 1 






s,t 



1 1 

S 
S 




1 

\ 2 
V 

\ 




dim s dim t 



+ s+t + 

_ s _t _ 



k 
(53) 

Again, to cover all possible orientation configurations, one simply needs to correct the formula 
by adding a phase (— l)^'' if a line of spin k is going in a direction opposite to the four others, 
and a phase (— 1)2('=+0 is the two lines k and I are pointing in different directions than the 
three others. 

Note that all the above integrals are invariant under permutation of the representations. 
Whilst a single intertwiner picks up a phase if two arguments are permuted, the whole integral 
is invariant because the phases appearing in each one of the two intertwiners appearing when 
performing the integrals cancel. 

Putting everything together, we obtain that the generating functional (|44p for a fixed 
assignment of spins factorises into amplitudes for the vertices of the dual triangulation 



<^"''^(j/,j)=n4"''^(j/'^)- 



(54) 



The vertex amplitudes are modified 6j symbols with source insertions |2] of the type 
discovered in [T]. Their explicit expression depends on the observable, and thus on the graph F. 
Typically, the vertex amplitudes are functions of the six spins and of the six sources assigned 
to the six wedges sharing the corresponding vertex. They are calculated by evaluating the 
representations on the sources. The indices of the resulting matrices tt-' (e*^) are contracted with 
the indices of the (three, four or five-valent) intertwiners issued from the Haar integrals. The 
tensor contraction pattern can be conveniently read out of the the diagrammatic representation 
of these amplitudes. 

The spin network diagrams associated to the vertex amplitudes are all built out of four 
vertices, or nodes, and eight (possibly nine) lines. The nodes represent the intertwiners issued 
from the Haar integrations and six out of the eight lines correspond to the spins coloring the 
faces, or wedges, meeting on the vertex v and contain a source insertion. The two (possible 
three) remaining lines are in the fundamental representation and are associated to the two 
edges of the graph F traversing the vertex or to one (possibly two) edge of F and one edge 
supporting the holonomy defining the observable if w is a vertex with a fermionic field 



^Noto that there could also be cases where a given edge e is associated to three lines, two issued from the 
graph r and one coming from a matrix Ve defining the gauge-invariant observable. In such cases, one uses the 
obvious generalisation of the above formulae. 



15 



insertion. Note that a matter line can form a loop based at a vertex if it is associated to an 
edge of r. In addition, the matter lines in the spinor representations are marked with a vertex 
with an open line in the adjoint representation colored by an adjoint index a or 6, weighted 
by the number i\/6/2, if they are associated to an edge of F. This is because to the edges 
e of r are assigned the quantities ja^Ue (see e.g. dH])). Thus the holonomies are composed 
with generators of su(2) in the spinor j = 1/2 representation. As explained below ([5^ . 
such a generator is an unnormalised intertwiner between the adjoint representation, a spinor 
representation and its dual and is thus represented diagrammatically by a vertex weighted by 
its norm, that is, the number iy/6/2. This closes our state-sum representation of the generating 
functional. We are now ready to calculate the action of the source derivative in ((M)) . 

Grasping operators. Deriving an exponentiated current in the j representation with re- 
spect to the a-th component of the source produces the evaluation of the su(2) generator ja 
in the Lie algebra representation j : 

= <{la), (55) 

.7=0 

where ni is the Lie algebra representation induced by the group representation tt^ . 

Now, the linear map iri : su(2) — > End Vj commutes with the group action and is therefore 
an intertwining operator between the adjoint representation, a j representation and its dual 

wi e Iloin{Vi,Vj(E)V*). (56) 

This interwiner is not normalised and its norm (squared) is easily evaluated to be 0(j)^ := 
<d{l,j,j)^ = + l)(2j + 1). Diagrammatically, a Lie algebra representation is thus 

represented as 




Therefore the action of source derivation with respect to the a-th component of the current, 
also called grasping, modifies a spin network with source insertion by creating a new vertex 
attached to an open line coloured by the index a in the adjoint representation. 

We can now readily see the action of the operator associated to an edge e of the graph 
r in p4p on the vertex amplitude associated to the vertex which is the source of e taken at 
J = 0. Each term in the sum over wedges in (j33| marks each one of the two representation 
lines associated to the two wedges We and w'^, with a vertex attached to an open line. The 
two resulting open lines are coloured by the indices b and c in the adjoint representation. 
These two lines, together with the open line associated to the edge e appearing in the vertex 
amplitude, connect through the totally antisymmetric three dimensional Levi-Civita tensor 
regarded as 3j symbol , 




+ 



This is the explicit picture of the soldering of the spin space to the quantum geometry of 
spacetime. One must then sum over over all admissible wedges. 

Finally, the amplitude ((M)) associated to a graph F and to the chosen configuration cr = 
is thus given by an expression of the type 

A" = {^a)^- n n E M'Jf.s.)- (59) 

Here, the sum over labels each of the elements of the bases of intertwiners for the intertwiners 
of valence greater than three appearing in the amplitude associated to the vertex v. A 
sum over graspings is implicitly assumed in the definition of the amplitude . Graphically, 
Ay is given by a (sum over) spin network diagram(s) weighted by a function of the spins 
N{i) obtained by assigning the quantity (— */4)0(j/e)0(i/^/ ) to each grasping occuring in 
the diagram, and a phase factor (—1)-^, with x the sum of all spins associated to the lines 
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oriented differently from the otlier lines at each vertex of the diagram which is issued from 
a Haar integration. All possible spin networks are given in Figure 1. Each one of these 
diagrams can be expressed in terms of 6j symbols using the recoupling identities given in the 
Appendix. For a particular triangulation and a particular graph/configuration, the obtained 
vertex amplitudes are recoupled in the following section. 

The extension to the other configurations c ^ associated to the graph F is immediately 
obtained from the above derivation. 

3.2.2 Gauge dependent observables 

When computing gauge dependent observables, the essential difference with the gauge invariant 
case comes from the fact that we can no longer suppress the gauge fixing function HeeT ^(ffe) 
from the computation if we do not want to obtain a vanishing expectation value. The resulting 
calculation will therefore depend on the gauge, and thus on the tree T. 

As in the gauge invariant case, we choose a particular configuration c — associated to an 
arbitrary graph F. The corresponding general form of the amplitude to be integrated yields 

4{e^„Ue)^{taf^Y[l[j:eUe lY[tiY[j:,uA. (60) 

V eeV \ C eeC ) 

All the steps performed in the gauge invariant case for the calculation of Aq transfer 
unchanged to the gauge dependent context up to equation (|44p. which is now modified with 
the inclusion of the gauge fixing function. After that, the derivation slightly differs from the 
gauge invariant case. The main consequence of the presence of the gauge fixing function is 
that not all group variables are integrated over. Therefore, although all formulae for the Haar 
integrals are implemented as in the gauge invariant case, the resulting generating functional 
at fixed spins (jH} can not put into factoriscd form as in 

Typically, (|54|) is replaced by a function of the spins and of the currents given by a collection 
of invariant tensors the indices of which are contracted with the indices of representation 
matrices evaluated on the currents. The pattern of contraction depends on the triangulation, 
the tree, the graph and can be potentially complicated. Thanks to the diagrammatic methods, 
one can read the contraction patterns out of the corresponding (open) spin network. The 
evaluation of this diagram is a function of the currents and yields the value of the generating 
functional 

<°''^(j/,^) = ((8)^^'(e-'0,(g)^n), (61) 

I n 

where I and n stand for the lines and nodes, or vertices, of the diagram, and the pairing (,) 
denotes the contraction of the indices determined by the combinatorics of the diagram. 

The implementation of the grasping operator proceeds in the same way than in the 
gauge invariant case and the amplitude associated to a graph F, to the chosen configuration 
cr = and to a particular tree T is thus given by an expression of the form 

< = (*«)^^ n E n E ^(^/' ^'')' (62) 

where ^(j/, St,) is given by the evaluation of a particular spin network diagram determined by 
the triangulation A, the tree T and the graph F. A sum over all possible graspings of each 
operator Se for each edge e is assumed in the notation so that A(jf, Sy) is in fact a sum over 
spin network diagrams. 

The exact form of A{jf, Sy) can be obtained by the following graphical method based on 
the cable and wire technology [28l [29] . 

• To each e of the graph F, associate an oriented open j = 1/2 line which is marked by a 
vertex from which emerges an open line in the adjoint j = 1 representation. 

• To each face / of the dual triangulation A* which does not touch any edge e of F associate 
a closed loop coloured by a spin j^j and orientated according to the orientation of the 
face. 
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Figure 1: Spin network diagrams for the different possible vertex amplitudes Ay (up to sym- 
metry) of the gauge invariant observables. The dashed lines are in the spin half representation 
and the curved lines denote the grasping operators in the spin one representation. To each 
depicted diagram, one must add the two diagrams obtained by summing over the possible 
grapsings. Note that due to the Ve terms there can be additional spin half lines that have not 
been included in these diagrams. We have also omitted for simplicity the ordering labels at 
each vertex and the arrows on the lines. 
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• For each edge e of the graph F, select a couple of faces /e, fe' (or of wedges We, We') among 
the three faces (wedges) touching e and s(e), and assign an oriented closed loop marked 
with a vertex from which emerges an open line coloured by an adjoint representation 
j = 1 to fe and fe' . 

• For each edge e of F, connect the two open lines emerging from the couple of faces fe, fe' 
to a three- valent vertex corresponding to the totally antisymmetric tensor 

• Connect, for each edge e of F, the third line emerging from the three-valent vertex to 
the open line emerging from the vertex marking the j — 1/2 line assigned to the edge e. 

• To each edge e of the complex A* not belonging to the tree T is assigned a box tra- 
versed by all the lines associated to e, which correspond to the irreducible representations 
labelling the wedges adjacent to the edge e possibly supplemented by spin 1/2 represen- 
tations 

• Perform the integrals by replacing the boxes by the corresponding intertwincrs using 
equations 

• Assign the factor (— i/4)0(j/J0(j/^, ) to each grasping 

• Sum over the diagrams obtained by summing over all possible couples fe, fe' for each 
edge e of F. 

The result of a such procedure for a given set of spins {jf, s^} yields a number A{jf, Sn) € C. 
From this spin network evaluation, one can calculate ■ The same prescription then applies, 
modified appropriately, to all configurations c of a given graph F. Note also that this graphical 
method applies to the gauge-invariant context by simply suppressing the tree. 

This closes the general prescription to compute observables in 3d spinfoam quantum gravity 
with fermions. To sharpen the details and techniques of the procedure, we now fix a particular 
triangulated 3-manifold and illustrate the techniques derived above to calculate the expectation 
value of the Polyakov line and of the fermion two-point function. 

4 Explicit calculations: observables on 

In this section, we give examples of the calculations described above for a space-time manifold 
homeomorphic to the three-sphere, i.e., M = S^. We consider the triangulation, denoted A, of 
S'^ with the boundary of a 4-simple3|3 and label each vertex v of A* (the pentachoron graph) 
by the labels I, J = 1, 5 which implies that the edges e of A* are labelled by couples IJ, 
with I ^ J, and the faces / of A* are labelled by triples UK, with I ^ J ^ K corresponding 
to the three vertices belonging to the face. The Dirac part of the action ^ for A, with rescaled 
triad, is 



The examples we will consider are the Polyakov line O fie) = (V's(e)j K;'/'t(e))i ^^nd the prop- 
agator O ^\(x\,X2) = ■0^1 i^x^B- We will now calculate the Berezin integrals explicitly to 
illustrate the Feynman rules for both observables. 

4.1 Gauge invariant observable: Polyakov line 

For this observable we consider 0/(23) = (i/)2, V253 V'3)i with V253 = V25V53. The bosonic 
observable O}, formed by integrating the Grassmann variables is 



®The order with which the hnes enter the boxes is irrelevant because of the permutation symmetry discussed 
above. 

''The (pseudo-)triangulation constructed from only two tetrahedra can not be used as it is not possible to 
properly define the action. This is because there is more than one edge connecting two vertices. While this 
is not problematic when working with theories of flat connections, it becomes highly ambiguous when dealing 
with non-trivial curvature. 




(63) 




(64) 
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We will consider one relevant non- vanishing term in the Feynman expansion (|19|) in order to 
compute the amplitude /r of a single admissible graph in equation (|21|) . 

X^1C-Di2dV'? V^2£-D23FV'f l^^icD^^H^" ^:iiDjM i^ARD^i^Li^i 
y-i^AMD^i^NlP^ ?5P^54qV'? V^5fl^51^S^f ?^1T^15 ^^s'' + (65) 

including the factor of 2 coming from the identical contribution of the iSfj and S^jj terms. This 
corresponds to a particular Feynman graph which we denote r(2, 3). Taking into account the 
order of the Grassmann variables, performing the integrals and contracting the indices leaves 
us with a global sign of er(2.3) = — 1- This agrees with the usual QFT Feynman rules where 
one associates a minus sign to each closed loop. 
Hence 

er(2,3)^r(2,3) = -(«a)^tr (£'32^253) tr {D21D15D54D41D15D54D43D32) ■ (66) 

Applying the graphical method described above gives us a number of admissible graphs, 
the graph r(2,3) corresponding to the term in the Feynman expansion above is given by the 
following 




The path 7^(3, 2) runs from vertex 3 to vertex 2 and has length E-p = 1, and the loop term C 
has Ec = 8 and passes through the following chain of vertices: 154154321. The overall length 
is therefore i?r(2.3) = 2n — p = 9. The Feynman rules in section [3] now associate a product of 
Dirac matrices De to the edges of r(2, 3) which gives equation 

This completes the fermionic integration, we can now consider the gravitational integration. 
As above, we will consider the c ~ configuration /p^j 3): which is given by 

Ir{2,3) = («a)'tr (£32(732^^253) 

xtr (S2l[/2lSl5?7l5S54C^54l]4iJ74iI]i5t/i5l]54C/54E43[/43) . (67) 

Using this, one can apply the graphical methods of section[3]to obtain the generating functional 
and then the amplitude Ap^2 3) of the graph r(2, 3). 

^r(2,3) = (f)' n 5]dimj,jK("l)'(^"=+^'^^^+^'^^=+-''^^^+-''^^=) (68) 

I<.I<KjlJK 

X E E E EEEEE 

Sl5,tl5 S23,t23 S45,t45 S14 S12 S25 S34 ^35 

X dim Sis dim ti^ dim S14 dim S12 dim S23 dim dim §25 dim S34 dim S35 dim S45 dim t45 

Xe2(j\23)e2(ji24)e(j\25)e2(j,,4)e2(j\35)e4(ji45)e(j234)e2(j235)e(j245)e(j345) 

X^l(jl.., Sl.)^2(j2.., S2.)^3(i3.., •S3.)^4(j4.., S4.)A(j5.., Ss.) 

Here, the global factor is due to the multiplication of the ia factors coming from the Feynman 
rules with the factor — i/4 coming from the global contribution attached to each grasping 
{iVQ/2 arising from the 7", a/6 from the e"'"^ tensor and —1/12 from the S"). 

The signs are due to the Haar integrals, and the 0{jijK) terms come from the normalisation 
of the source derivatives. The vertex amplitudes Aj, I = 1, 5, are given by evaluating the 
following spin networks. 



20 




21 



start by decomposing the vertices of valence four and five into tlirce-valent vertices. Then, we 
repeatedly use the recoupling identities ([98|) . (jlOOp . (|102p displayed in the Appendix. As a 
result, we obtain the following decomposition of the vertex amplitudes Aj{jjjK, s/j) appearing 
in equation 




A. 




^ (dim S23 )" Mini m dimr (5^23,525 (-1)^^^23+3^245+^234) 



X(-l)("^+2^ + 2^Si2+2ji23+j235+3i23+3il25+Jl24 + l) 



Jl25 Jl24 ^23 

J234 -^23 J245 

Jl25 ^ Jl24 

1 ii24 r 






, J234 




^23 , 




■J345 \ 


^3 = S23~. 




J123 / 




■"C^*. .7235 

— TiN^ 


/0'l35 



^^p^^dimri dimp dimg dimu I)(ii35+j345+s34) 

X l)(3" + 9+"+jl23 +2^234 + 2*23 +S23+il23+S35 +2 jl34+S23) 



A. 




^^j^dima dim 6 dime 

X ( — I)(3a+3b+3c+3s34+j234+jl45+j345 + f +3S45- 



-3jl34) 



.7l34 J234 ^45 

,7245 ^45 J124 

1 1 1 

a c 

-^14 J134 J124 
Jl24 1 & 



J134 J234 ^45 
C 1 jl34 

1 1 1 

14 & 

1 ^14 ^4 
7l45 2 2 



S45 



J134 J124 -^45 
1 jl34 
^ 7l34 J234 

1 J234 c 
J234 2 

J345 J 134 
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-^.7145 \ 




J345 \. 











+ 
S3 
+ 



Ea,b,c,d,ej,g,h,»dimadimfe dime dimd dime 
xdim/ dim(^ dim/i dimi ^S23,s25(dims3g)~^ 

x(-l)(3c+2d+2/+23+3/i+3i+12ji2i+3i245+3s25) 

X ^_ 2^ ^ (3 jl45 +il35 +3i345 +3si5 +3^45 +i235 +S35 ) 

J125 J245 'S25 11 * J345 J24 11 1 j 
J345 jl35 i J \ j345 3 1 J I .?145 ^ 
J125 ^ 1 / 3 \ / ^ ^ ^ 

ji45 ji35 « J \ ji45 1 e J \ e d /l 

a e 1 1 f 1 1 
J145 ji45 g ] \ a e 



! 
2 
h 



1 e 
]_ ]_ 

2 2 

J345 ^4 
]_ 
2 



a 



2 

^45 
\_ 
2 



Jl4 

e 

1 

1 

2 '^15 
2 d 

J245 ^4; 

jl45 ^ 



1 a 

1 1 

2 2 

J135 -^IS ^1 

2 J145 C 



This completes the evaluation of 3)- 

4.2 Gauge variant observable: Fermion propagator 

For the propagator we calculate one term of the expectation value of C'y"^(2,3) = 1^2 

in a particular gauge. For this observable, the term in the Feynman expansion (|19p that we 

consider is 



Ob 



+ {iaf 
xV'ic-Di2dV'? V' 



dfl^Tpj, tpl)] ■02 ^. 



^33 



2£^21FV'f VlG^ls'^fV'.f ^S/^SlVV'l' V' 
34ArV'4'' V'4P-D43qV'? ^/'4fl£'45^V'f V'5T^54^V'4 



/,C/ 



2K^23^'/'3 



(74) 



including the factor of 2 coming from the identical contribution of the and S^jj terms. 

Applying the graphical method gives the admissible graph r'(2,3) corresponding to the 
above term in the Feynman expansion 





This gives a single path V'{2, 3) of length i?r'(2,3) = E- 



•P'(2,3) = 9 passing through the vertices 
2154321543. Performing the integration and contracting the indices as before gives the sign 
factor er(2,3) — 1 since there is no closed loop. The amplitude /r'(2.3) of the graph r'(2,3) is 
given by 

^r'(2,3f B = {D2iDi^DziDM2D2iDi^D^iDiz)% . (75) 

This completes the fermionic integration, we can now consider the gravitational integration. 
As above, we will consider the c = configuration /p,(2 3)1 which is obtained from the above 
expression by replacing the Dirac matrices Dij with the quantities YiijUij. 

We then choose the maximal gauge fixing tree T on A* given by the set of edges T = 
{21,15,54,43}. 

We can now apply the graphical methods of section |3] to obtain the generating functional 
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and then the amphtude A^,^^ 3) of the graph r'(2,3) for the tree T. We obtain 

^r'(2,3) = (f)^ n ^ dinii/j/f(-l)2(Ji«+Ji2=s+J235+j234+ji34)+i (7g) 

K.KKjijK 
X diniS23 

S23 

X o\,3 ) (ji25 ) (.7134 ) (j,,, ) 0,34 ) O345 ) 

where the amphtude jat, S23) G C is obtained by evaluating the spin network diagram 
displayed below. 



A{jlJK,S23) = 




By Schur's Lemma, this will be proportional to the identity map Sg in the fundamental 
representation multiplied by half the value of the closed spin network evaluation obtained by 
connecting the two open lines in the j = 1/2 representation. The closed diagram can then be 
recoupled and given as a combination of 6j symbols. 

5 Conclusions 

In this paper, we have shown how to compute expectation values of fermionic observables in 
spinfoam models of 3d quantum gravity. We have identified and properly gauge fixed the 
symmetries of the model, given a general procedure to compute gauge invariant and gauge 
dependent observables and illustrated the techniques on an example. 

It is important to stress the achievements of this work with respect to the original proposal 
[1]. We have 1) extended the model to massless fermionic fields, 2) implemented a new 
procedure to compute the Berezin integrals, 3) introduced observables, 4) recoupled the vertex 
amplitudes to express them as combinations of 6j symbols. 

Many issues remain open. First, unlike the point particle case, fields break the topological 
invariance of three dimensional gravity. The model is therefore only an approximation of the 
continuum theory, even if the spinfoam explicitly derives from a classical action. This leads 
to two questions. First, do the asymptotics {j 00) of the vertex amplitudes reproduce (a 
simple function of) the classical Regge action coupled to the Feynman diagrams issued from 
the fermionic integration? Second, and most importantly: what is the behaviour of the model 
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under refinement of the triangulation? Since the finiteness of our Fcynman diagram expansion 
is due to the cut-off on the number of degrees of freedom induced by the triangulation, this 
issue is of primary interest. This issue was considered in |30| where a group field theory 
approach to the model was constructed, in order to implement a sum over triangulations. 
Finally, one would need to analyse the well known fermion doubling problem appearing in 
lattice-like theories incorporating fermions. Does the spinfoam generate doublers? Could we 
modify the action to suppress the extra species? 

We close with the perspectives. Working with a three-dimensional non-topological model 
has taught us techniques that could be useful in the full four-dimensional theory coupled to 
matter. This leads us to the most interesting outlook which would be a generalisation of the 
ideas and techniques developed here to four spacetime dimensions. At first, one might expect 
the coupling of fermions to the EPRL or FK models [311 E21 El] of four dimensional quantum 
gravity to be problematic. The reason is that the kinetic term of the standard Dirac action in 
a d-dimensional curved spacetime couples to the d — 1 exterior power of the d-bein. Hence, in 
even dimensions, fermions are coupled to Palatini gravity via an odd power of the gravitational 
field not convertible into any power of the d — 2-form B field of BF theory. However, if one 
uses the Capovilla, Dell, Jacobson, Mason chiral action [34] of fermions coupled to two-form 
gravity, the fermions indeed couple to the B field and one could imagine such an action as an 
interesting starting point for a four-dimensional construction. 
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A Basic facts on the representation theory of SU(2) 
A.l Notation and conventions 

Our conventions are essentially those of references [35l |36]. The unitary irreducible repre- 
sentations of SU(2), ttJ : SU(2) Aut(Vj), are labelled by a spin j e N/2. The spin j 
representation acts in the finite dimensional complex vector space Vj = C^-'+^. In fact, Vj can 
be endowed with a metric structure inherited from the pairing map : Vj —t' V* , where V* 
is the vector space of linear forms on Vj. Choosing a particular basis {el^)a=-j....,j of Vj, the 
components of the paring map, the Wigner I j symbol, are given by 

This map satisfies the properties 

The associated inner product is defined by {v,w)j = e^{v){w) = v^w^e^^p := v^Wa, for v^w in 

For all triple of unitary, irreducible representations tt', tt-', tt*^ of SU(2), let 

C% : V ® V, ^ Vk, and C% : Vk ^ V, (S) V,-, (78) 
denote the Clebsh-Gordan intertwining operators whose coefficients 



: ^ ) ei;, (79) 



are well known from the quantum mechanics of angular momentum. These intertwiners are 
unique up to normalisation because the space of three-valent intertwiners is a vector space of 
dimension one. The complete reducibility of the tensor product of representations 



V^V,^ Vk, 

k=\i-j\ 
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is then expressed as tt* (g) tt^ 



J2k ^% ^'^ Relative to the ehosen bases, this yields 



E 

fc,e,C 



a 7 



P 6 



(80) 



The Clebsh-Gordan coefficients satisfy rather awkward symmetry properties . One can 
consider a rescaled version of these symbols that are more symmetric. These objects are 
called 3j symbols. They are evaluations of the 3j maps which are normalised three-valent 
intertwiners. Let Lijk :Vi®Vj®Vk ^ 'Che & 3j intertwining map and w, v and w be arbitrary 
vectors in Vi , 14 and 14 respectively. Clebsh-Gordan symbols and 3 j symbols are related via 
the following evaluation 



iijk{u ®v®w) 



V dim k 



(81) 



where e(«, j, k) is a sign usually chosen to be (—1) ^ ^ where —i+j~k G Z. The evaluation 
of the above equation on the vectors u — e'^, v - 



and w 



reads 



i j k 

a j3 ^ 



Vdimfc 



a P 



(_l)-»+J+7 / _^ 

%/ dim fc V 



« .7 



(82) 



where the left hand side is a 3j symbol. The above relation can be put into a more compact 
form by using the symplectic metric to raise and lower indices. Our conventions are such that 
indices are raised from the right and lowered from the left, that is, picturing only one column 
of the symbol, 



/3 



and 



Note that this implies the following property 



(-1) 



2i 



A Pa 



(83) 



From these rules, one can construct mixed 3j symbols, that is, 3j symbols associated to all 
possible three-fold tensor products involving three vector spaces together with their duals. 
Furthermore, using the symmetry under permutations of the 3j symbol 



i j k 

a /3 7 



(-1) 



i k j 

a "f 13 



(84) 



which implies that the symbol is even under even permutations of the columns, the relation 
between 3j symbols and Clcbsch-Gordan coefficients can be put in the form 



k i j 

J a (3 



V dim k 



« J 

a (3 



(85) 



A. 2 Group integrals 

Integrals of tensor product of representations can be expressed in terms of 3j symbols. The 
first step is to realise that 

SU(2) 

From the above expression and the complete reducibility, one can compute all integrals involved 
in this paper. The first in the bivalent integration. The calculation proceeds as follows. 



SU(2) 



7 _ 



E 

fc,e,C 



a 7 



^ J 
P 5 



dgTT^{g) 



c 



SU(2) 



/ a 7 


S) 


(o 


i 





V * 3 






/3 





iim J 



(86) 
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where we have used that 



a p 
i j 



\/dim j 



in the last step. From this resuh it is immediate to prove the orthogonahty of characters ()4ip . 

Using the above, one can compute the three-valent integral. The calculation is displayed 
below. 



SU(2) i,-,^g 



a 7 
i j 



I 



1 



dim k 



ri,e 



a 7 
i j 



« J 

/3 S 



keek 



SU(2) 



dgiT\g)\','^{g)\ 



/3 5 



/ a 7 e \ 

















(87) 



To calculate the same integral, but with one group element inverted, on simply uses the 
fact that the dual pairing is in fact an isomorphism of representations, that is, a bijcctive 
intertwiner. This translates as 



From the above, we immediately obtain that 



(88) 



dg^\grpn={g)\^\g'')\ = {^ir 



\2k 



SU(2) 



a 7 A: 
i J C A /3 5 k 



I J e 



Because of the invariance of the Haar measure under inversion g i— > g^^, we have, with the 
two above formulae, covered all possible cases of orientations that can occur in the model. 
We can next proceed analogously for the four-valent case. We obtain the following result 

/ dg7:\gr^7:^igyy{gy^n\g)%= ^ ^ 



/ Of 7 


m \ 
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V * j 
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\ m 
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dim™ 
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E 



a 7 e 
i j k I 
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i j 



k I 
C V 
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SU(2) 
m i.^ 



dg7r"^{gY^n-{g)\ 



k I 
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« J 
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k I 



i j k I 

P S C Tl 



(89) 



where the obtained 4j symbol is defined as follows 
a P 1 5 ) ^ 



I J K 

a P m 



m k I 
K ^ 6 



(90) 



Using the conjugation (|88p . we can also calculate the following integral 
dgn^igr,n^{gr,7:\gr,.\g~r, = i-irY.('t ] I \ 



/SU(2) 

and also 



i j k 9 
(3 S C I 



SU(2) 



rf.9-'(5)",-^(ff)V'=(5-^)\-'(.g"^)% = (-i)^('=+')E(" ] c I 



I ] e 
fi 5 k I 



This covers all possible four-valent orientations configurations. 
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Finally, we compute the five-valent integration as follows. 
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(91) 



where the 5j symbol is given by the following expression 

i j jjL \ I n k a 



i j k I m 
(3 5 Q T] K 



= \/ dim n dim ( 



13 5 n J \ jjL ( o 
Again, from the above, it is immediate to obtain that 

dgi:\grf,i:\g)\i^\g)\,^\g)\i:^\g-^)\ 
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and that 



SU(2) 



dg^\grpir={g)\^\g)\^\g-^)\^^{g-^)\ 
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a 7 e Z m 



i j k 9 L 
13 5 C I m 



i j k rj K 

This covers all possible group integrals that can appear in the model. 



A. 3 Diagrammatics and Recoupling identities 

The spin network vertex amplitudes obtained in this paper can be re-expressed as products 
of 6j-symbols using recoupling identities. These identities are most conveniently express dia- 
grammatically. A 3j symbol is represented as an ordered vertex with three lines. The clockwise 
(resp. anti-clockwise) ordering of the lines at a vertex is denoted with a -I- (resp. — ) symbol. 
If required, arrows will be depicted on the lines with the convention that ingoing lines at a 
vertex correspond to lowered magnetic components in the symbol. This implies that 



i j k 

a (3 J 



and 



a /3 7 
i j k 



The 3j symbols are normalised 



i j k 

a /3 7 




a /? 7 
i j k 



Y{i,j,k), 



(93) 
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where Y{i^j,k) = 1 if i,fc and / satisfy the triangle inequahties and yields zero otherwise. 
Diagrammatically, this means that the theta network is equal to (zero or) one 





(94) 



k 



Note that changing the orientation of an arrow, say of spin k, in the above diagrams produces 
a factor (— l)^*"'. The corresponding intertwiner has norm (squared) given by (~l)^'^K(i, j, k), 
in accordance with the results found in the group integral section. 

Changing the ordering at a vertex is equivalent to a permutation of the arguments of the 
3j thus we have the diagrammatic identity 




\« J / 



(95) 



This emphasises the need to have an ordering label at the vertices. 

The 6j symbol is defined by the following four-fold contraction of 3j symbols 
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ji k je 
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Accordingly, the 6j symbol is equal to the following tetrahedral spin network evaluation 



31 32 33 

34 J5 36 




(96) 
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Note that reversing an arrow on a spin j line of the diagram produces a factor (—1)^-' because 
of the 'see saw' property ([55)1 of the inner product. This emphasises the need to have the 
arrows in the diagram. Thanks to the symmetries of the 3j symbol, the 6j symbol network is 
invariant under any permutation of the vertices. 

We are now ready to introduce the recoupling properties used throughout the paper. The 
first that we use is the second orthogonality of the 3j symbols 



i 7 (5 
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k I /3 
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5) 5iY{i,k,l). 



(97) 



This property can be understood as Schur's lemma. This is most transparent in the diagram- 
matic version of the above identity 
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(98) 



We also made use of the recoupling theorem which provides the change of basis formula 
between two different recoupling schemes of a 4-valent intertwiner 
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where 



Cn = dim 71 (— 1 



yi+k+m+n J « J 771 

k I n 



29 



In graphical language, the recoupling theorem yields 




A consequence of the recoupling identity is the following property 

f i 6 n \ f I i e \ ( m k ^ \ _ i i k j \ ( i j k \ , . 

\ a I e J \ S P m J \ 7 n J \ m I n j \ a P 77" ^ ' 

Diagrammatically, this yields the fusion move displayed below 
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